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Abstract 

A  general  paradigm  for  recognizing  3D  objects  is  offered,  and  applied  to  some  geometric  primitives 
(spheres,  cylinders,  cones,  and  tori).  The  assumption  is  that  a  curve  on  the  surface,  or  a  pair 
of  intersecting  curves,  has  been  measured  with  high  accuracy  (for  instance,  by  a  sensory  robot). 
Differential  invariants  of  the  curve(s)  are  then  used  to  recognize  the  surface.  The  motivation  is 
twofold:  the  output  of  some  devices  is  not  surface  range  data,  but  such  curves.  Also,  a  considerable 
speedup  is  obtained  by  using  curve  data,  as  opposed  to  surface  data  which  usually  contains  a  much 
higher  number  of  points. 

We  survey  global,  algebraic  methods  for  recognizing  surfaces,  and  point  out  their  limitations. 
After  introducing  some  notions  from  differential  geometry  and  elimination  theory,  the  differential 
and  “semi-differential”  approach  to  the  problem  is  described,  and  novel  invariants  which  are  based 
on  the  curve’s  curvature  and  torsion  are  derived. 
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1  Introduction  and  Previous  Work 


One  task  an  intelligent  system  should  be  able  to  accomplish  is  recognition.  Usually,  a  recognition 
system  derives  some  characteristics  of  an  object  it  examines,  and  tries  to  match  them  against  similar 
characteristics  in  a  data  base.  Suppose,  for  instance,  that  one  is  dealing  with  2D  objects,  and  tries 
to  recognize  them,  given  their  boundaries.  Typically,  there  is  a  finite  data  base  these  boundaries 
are  matched  against;  various  invariants  have  been  derived,  some  global  and  some  local  [1,  17,  26], 
to  solve  this  problem.  These  are  quantities  that  do  not  change  under  certain  transformations 
(Euclidean,  affine,  projective),  and  therefore  can  be  used  to  recognize  an  object  even  if  it  has  been 
altered  by  such  transformations. 

Here,  a  different  problem  is  addressed  -  recognizing  a  surface  in  3D  space,  using  one-dimensional 
information.  Specifically,  we  assume  that  some  measuring  device  has  sampled  a  curve,  or  a  pair  of 
intersecting  curves,  on  the  surface.  Given  the  curve(s),  the  goal  is  to  recognize  the  surface.  Typical 
sensors  which  are  the  source  of  such  curves  are  measuring  devices,  such  as  coordinate  measuring 
machines,  manufactured  by  the  Brown  k  Sharpe  Company  (Figure  1),  or  the  IBM  RS/l  Cartesian 
robot.  Such  devices  can  measure  3D  curves  with  very  high  accuracy;  for  instance,  the  typical  error 
range  for  a  coordinate  measuring  machine  is  0.01  mm. 

In  [3],  an  algorithm  is  presented  for  determining  the  axis  of  a  surface  of  revolution,  using  the 
information  measured  by  a  tactile  sensor  which  can  also  estimate  the  two  principal  curvatures  (see 
Section  3.2).  In  [7],  the  parameters  of  a  cylinder  are  computed  from  structured  light  patterns. 

Some  previous  work  has  addressed  the  problem  of  recognizing  various  surfaces  given  their  oc¬ 
cluding  contours  [14,  10].  However,  the  aggregate  of  possible  curves  on  a  surface  is  usually  much 
larger  than  the  aggregate  of  its  occluding  contours,  and  may  contain  far  more  complicated  curves; 
for  instance,  the  occluding  contour  of  a  sphere  is  always  a  circle,  while  there  are  a  great  many  3D 
curves  -  some  of  which  have  rather  complicated  structures  -  on  a  sphere. 

Clearly,  we  are  facing  a  different  type  of  recognition  problem  from  the  one  previously  described, 
which  is  usually  solved  by  matching  against  a  data  base.  It  is  impossible  to  build  a  data  base  which 
contains,  say,  all  the  curves  on  a  sphere,  or  even  a  dense  sampling  of  these  curves.  Therefore,  we 
have  to  discover  curve  characteristics  which  will  enable  us  to  answer  a  question  such  as  “can  this 
curve,  after  a  certain  transformation,  be  embedded  into  a  sphere?”,  as  opposed  to  “can  this  curve, 
after  a  certain  transformation,  be  superimposed  onto  curve  No.  129  in  the  data  base?”. 

One  way  to  proceed  is  straightforward:  fit  an  implicit  polynomial  to  the  curve’s  points,  and, 
from  its  type,  determine  the  surface.  This  is  the  algebraic  approach  [19, 12].  However,  this  approach 
win  fail  if  the  curve  does  not  lie  on  a  single  “primitive”  (sphere,  cylinder  etc),  but  “crosses  over” 
between  two  or  more  primitives  (see  Section  2.1  and  Figure  3).  In  that  case,  the  global  algebraic 
fit  will  give  us  a  meaningless  result.  A  very  rich  theory  of  local,  or  differential,  invariants  was 
developed  to  solve  this  problem  [5,  4,  8,  24,  21].  In  Section  2  we  quickly  survey  the  global  approach 
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Figure  1:  High- accuracy  measuring  device  and  a  curve  it  measured  on  a  cylinder. 


as  applied  to  our  problems,  but  the  focus  of  this  paper  is  on  the  local  approach. 

Natural  curve  characteristics  to  use  for  recognition  are  curvature  and  torsion,  as  they  do  not 
change  under  rigid  transformations.  Since  we  are  dealing  with  3D  data,  a  rigid  transformation 
is  usually  a  general  enough  model.  So,  the  goal  is  to  discover  invariants,  depending  on  a  curve’s 
curvature  and  torsion,  which  provide  a  necessary  condition  for  it  to  lie  on  a  certain  type  of  surface. 

Let  us  demonstrate  this  by  a  simple  2D  example:  a  plane  curve  can  be  embedded  into  a  circle 
if  and  only  if  its  curvature  is  constant.  So,  in  this  case,  the  invariant  is  the  derivative  of  the 
curvature.  We  don’t  expect  to  find  such  simple  invariants  for  curves  lying  on  3D  surfaces;  a  well- 
known  example  is  the  fact  that  a  3D  curve  is  planar  if  its  torsion  is  zero,  but  this  is  an  exceptional 
case. 

In  the  sequel,  we  derive  invariants  which  are  necessary  conditions  for  a  curve,  or  an  intersection 
of  two  curves,  to  lie  on  a  sphere,  cylinder,  cone,  or  torus.  These  depend  only  on  the  curvature 
and  torsion  at  a  point  on  the  curve  (or  the  curvatures  and  torsions  of  the  two  curves  at  their 
intersection  point).  We  also  derive  some  “semi- differential”  invariants,  which  use  not  only  the 
differential  properties  of  the  curve,  but  a  few  points  on  it.  Such  invariants  have  been  widely  used 
in  computer  vision  for  recognizing  plane  and  space  curves  [4,  21,  18];  their  main  advantage  is  that 
they  allow  us  to  use  derivatives  of  lower  order  than  are  needed  for  “purely  differential”  invariants. 
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2  The  Algebraic  Approach 

Implicit  polynomials  can  be  used  to  describe  2D  and  3D  objects.  Some  references  which  address 
the  fitting  of  implicit  polynomials  are  [22,  2,  13,  23,  20].  One  can  then  use  polynomial  invariants 
to  recognize  the  objects  [19,  12,  9,  11].  Let  us  briefly  describe  how  a  sphere,  cone,  cylinder  and 
torus  can  be  recognized  using  such  invariants.  Note  that  the  first  three  objects  can  be  fitted  with  a 
quadratic,  and  the  torus  with  a  quartic.  Suppose,  then,  that  we  have  succeeded  in  fitting  the  data 
with  a  quadratic.  We  can  write  this  as 

+  (v,  X)  +  s  =  0  (1) 

where  A  is  a  3  x  3  matrix,  v  a  vector  in  and  s  a  scalar.  It  is  easy  to  verify  that 

•  If  the  object  is  a  sphere,  A  has  three  positive  and  identical  eigenvalues.  It  is  then  trivial  to 
extract  the  sphere’s  center  and  radius. 

•  If  the  object  is  a  cylinder,  A  has  two  positive  and  identical  eigenvalues,  and  one  zero  eigen¬ 
value;  also,  the  axis  of  the  cylinder  is  in  the  direction  of  the  eigenvector  with  zero  eigenvalue, 
and  it  is  trivial  to  extract  its  radius. 

•  K  the  object  is  a  cone,  A  has  two  identical  positive  eigenvalues  and  one  negative  eigenvalue. 
The  axis  of  the  cone  is  in  the  direction  of  the  eigenvector  with  the  negative  eigenvalue.  It  is 
then  trivial  to  extract  the  cone’s  opening  angle  and  apex. 

•  If  the  object  is  a  torus,  its  general  equation  is 

Etor  =  {{x  -  af  +  {y-  hf  +  (z  -  cf  ^  -  4  R\{x  -  af  + 

{y  -  bf  +  iz-  cf  -  ((x  -  a)ni  -h  (y  -  h)n2  +  (2  -  c)nsf) 

where  (a,6,c)  is  its  center  point,  (ni,  712,^3)  is  a  unit  vector  perpendicular  to  the  plane  over 
which  the  torus  lies,  and  R  (r)  is  the  major  (minor)  radius. 

It  is  trivial  to  extract  a,  6,  c  from  Etor  (for  instance,  differentiating  Efor  three  times  wrt  x 
gives  24a:  -  24a).  To  extract  r  and  J2,  note  that  substituting  {x  =  a^y  =  b^z  —  c}  in  Etor 
gives  -1-  -  2  and  substituting  {x  =  a,  y  =  b,z  =  c}  in  ^ 

— 12  —  12  -|-  8  -|-  8  R^n2^  -|-  8  R^nf  =  —4R^  —  12r^.  It  is  trivial  to  extract  R  and 

r  from  these  two  identities.  After  R,  r,  a,  6,  c  have  been  recovered,  it  is  trivial  to  recover 
(711,712,^3). 


2.1  Number  of  Points  Needed 

Experiments  on  curve  data  show  that  a  relatively  high  number  of  points  is  necessary  to  achieve 
reliable  algebraic  fitting.  For  instance,  for  the  cylinder  data  we  have  used  (Figure  1),  more  than 
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200  points  are  required  for  a  reliable  fit.  We  are  not  sure  why  this  happens;  apparently,  the  fact 
that  the  points  lie  on  a  curve,  which  is  a  “one- dimensional  entity”,  results  in  singularities  when 
trying  to  fit  the  points  with  an  implicit  polynomial  which,  by  its  nature,  is  appropriate  for  fitting 
“two-dimensional  entities” . 

On  the  other  hand,  when  using  the  differential  invariants  proposed  here,  a  far  smaller  number 
of  points  was  necessary;  usually,  invariants  were  computed  using  10  points  or  so. 


2.2  Applying  Invariants  to  Segmentation 


Since  the  algebraic  approach  to  recognition  given  a  curve  may  fail,  because  it  can  pass  through 
various  geometric  primitives,  one  may  try  to  segment  the  curve,  using  some  notion  of  discontinuity, 
and  then  use  algebraic  techniques  for  each  segment.  We  now  show  that  this  is  not  always  easy,  by 
constructing  a  curve  which  is  infinitely  differentiable,  yet  crosses  over  from  a  sphere  to  a  cylinder. 
Define 


s{t)  = 


0  t<0 

exp('^)  t  >  0 


It  is  weU  known  that  this  function  is  smooth  (infinitely  differentiable)  at  every  point,  and  that 
all  its  derivatives  at  t  =  0  are  zero.  Using  s(t),  it  is  trivial  to  construct  smooth  functions  Si{t),S2{t) 
on  the  interval  [0,  oo)  such  that  Si(0)  =  0,  Si{t)  =  1  for  t  >  1,  S2{t)  =  y/Z  for  0  <  t  <  1,  and  S2{)  is 
monotonically  increeising  for  t  >  1  (see  Figure  2). 


Figure  2:  Auxiliary  functions  used  to  construct  the  curve  in  Figure  3. 
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Define  a  curve  c(i)  as  follows: 


ct  = 


{si{t)  cost,si{t)  sint.yjA  -  s\t  )  0  <  ^  <  1 
(cos sin S2{t))  1  <t  <2 


It  is  easy  to  see  that  c(t)  is  a  smooth  curve  which  crosses  over  from  a  sphere  with  radius  2  to 
a  cylinder  with  radius  1  (at  ^  =  1).  The  curve  is  displayed  in  Figure  3.  Next  to  it,  we  plot  the 
curvature,  torsion,  derivative  of  the  curvature,  and  a  spherical  invariant  for  curves  (see  Section  5, 
Equation  10).  It  is  interesting  to  see  that,  although  the  curvature  and  torsion  are  continuous,  there 
is  a  very  sharp  break  in  the  spherical  invariant,  at  the  point  at  which  the  curve  crosses  over  from 
the  sphere  to  the  cylinder;  this  demonstrates  that  the  kind  of  invariants  presented  here  can  succeed 
where  segmentation  by  ‘^ordinary”  differential  properties  (curvature,  torsion  etc.)  fails. 


0.2  0.4  0.6  0.8  1  1.2  1.4 

Figure  3.  Demonstrating  how  our  invariants  manage  to  detect  when  a  curve  crosses  over  from  one 
geometric  primitive  (sphere)  to  another  (cylinder),  although  the  curvature,  torsion  etc.  cannot 
detect  this  crossing  over. 
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3  Mathematical  Preliminaries  -  Some  DilFerential  Geometry  and  Elimination  Theory 


In  the  sequel,  a  few  concepts  from  geometry  and  algebra  are  required.  We  proceed  to  define  them 
and  state  some  of  their  important  properties. 

3.1  Some  Differential  Geometry  of  Curves 

A  curve  in  3D  Euclidean  space  is  a  differentiable  function  c  :  [0, 1]  -+  TZ^.  At  each  point  c(t),  three 
orthogonal  unit  vectors  are  associated  with  the  curve:  its  tangent  vector  T,  which  points  in  the 
direction  of  the  curve’s  derivative,  its  normal  vector  N,  and  its  binormal  vector  B,  which  is  equal 
to  the  vector  (cross)  product  of  T  and  N.  This  triplet  of  vectors  is  called  the  Frenet  trihedron  at 
c(t). 

In  addition,  two  scalars  are  associated  with  each  point  on  the  curve.  These  are  the  curvature  k 
and  torsion  r.  Intuitively  speaking,  the  curvature  measures  how  ^^bent”  the  curve  is;  for  instance, 
the  curvature  of  a  circle  is  equal  to  the  inverse  of  its  radius.  The  torsion  measures  the  speed  at 
which  the  curve  moves  out  of  the  plane  (the  so-called  osculating  plane)  which  locally  approximates 
it;  thus,  the  torsion  of  a  planar  curve  is  zero. 

The  curvature  and  torsion  can  be  computed  from  the  parameterization  of  the  curve: 

^  _  |c  X  C  I 

T 

K  and  r  are  invariant  to  translation  and  rotation;  this  makes  them  especially  attractive  for  recog¬ 
nition  purposes. 

The  celebrated  Frenet  formulas  relate  the  Frenet  trihedron  to  the  curvature  and  torsion.  If  the 
curve  is  parameterized  by  arclength  (that  is,  |c  |  =  1),  the  following  hold: 

N'  =  -kT  -  tB 


B'  =  tN 


A  concept  of  crucial  importance  to  our  work  is  the  local  canonical  form.  Let  us  see  how  it  is 
derived.  Assume  that  the  curve  is  parameterized  by  its  arclength  5.  From  Taylor’s  expansion,  we 
have 


c(«)  =  c(0)  +  sc'(O)  +  |ic"(0)  +  ^c'"(0)  +  ^cM)(0)  +  |;c(»)(0)  +  <,(**) 
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c'(0)  is  equal  to  the  tangent  vector  T  at  t  =  0.  Using  the  first  Frenet  formula,  c"(0)  =  T'  =  kN. 
Therefore,  c"'(0)  =  (kN)'  =  k  N  +  kN'  =  kN  +  k{-kT  -  tB)  =  kN  -  k^T  -  ktB. 

Similarly,  we  can  derive  expressions  for  the  fourth  and  fifth  derivatives.  Substituting  them  into 
the  Taylor  series  gives 


t  \  I  -T-  I  N  -  T  -  KT  B) 

c{s)  =  c(0)  +  sT+  - - - 1  + 

s^(k"  N  -  2k  t  B  -Zkk  T  -  k^N  -  kt  B  -  kt'^N) 

24 

S®  m 

N  -  a  kk"  T-ZktB-Zkt'B-Zk  t^N  -Zk^T  - 
6 K^K  N  k^T  A-  K^T B  -  kt"  B  -Zktt  N  k^t^ T  +  kt^B)  +  o(.s®) 


(2) 


From  now  on,  we  shall  omit  the  0(5®)  part.  We  are  allowed  to  do  so  as  long  as  the  powers  of  s  used 
are  bounded  by  5. 


3.2  Some  DiiFerential  Geometry  of  Surfaces 

Locally,  a  surface  5  in  3D  Euclidean  space  is  a  differentiable  image  of  an  open  set  in  7^^. 
Formally,  it  is  the  set  of  triplets  {{x(^u,v),y{u,v),z(^u,v))  (u,  u)  €  O}.  The  tangent  plane  to  S  at 
the  point  ((x(u,  v),  y{u,  u),  z(u,  u))  is  the  plane  spanned  by  (x^,  yui  ^u)  and  (x^,,  y^,  z^).  The  normal 
to  S  at  (u,  v)  is  the  unit  vector  pointing  in  the  direction  of  (x^,  y^^,  z^)  x  {xy,  y^,  z^);  it  is  therefore 
perpendicular  to  the  tangent  plane. 

In  the  sequel,  we  shall  use  the  fact  that  if  Oi  and  O2  are  curves  which  intersect  on  S,  then  the 
normal  to  5  at  their  intersection  point  is  a  unit  vector  in  the  direction  of  the  vector  product  of 
their  tangent  vectors.  This  holds  unless  these  tangent  vectors  are  parallel. 

The  intersection  of  S  with  any  plane  containing  N  is  called  a  normal  section  of  S.  Note  that 
the  normal  section  is  determined  by  a  unit  vector  v  in  the  tangent  plane,  which  is  the  direction 
in  which  the  plane  containing  N  intersects  the  tangent  plane.  Thus,  we  may  speak  of  a  normal 
section  in  the  direction  v. 

The  curvature  of  a  normal  section  is  called  the  normal  curvature.  The  maximal  such  curvature, 
ki,  and  the  minimal,  ^2,  are  called  the  principal  curvatures  of  S.  Let  us  denote  their  directions  by 
ki  and  k2.  It  can  be  proved  that  they  are  orthogonal  and  that,  if  u  =  cos  0  sin  0,  the  normal 
curvature  in  the  direction  v  equals 

ki  cos^  0  +  k2  sin^  0  (3) 

The  product  R  =  k\k2  is  called  the  Gaussian  curvature.,  and  the  mean  H  =  is  called  the 

mean  curvature. 

Suppose  a  curve  C  lies  on  the  surface  S.  Then  if  its  curvature  is  kq,  and  the  normal  curvature 
of  S  in  the  direction  of  C’s  tangent  vector  is  K5,  we  have 

xs  =  xc  cos  6 
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(4) 


where  9  is  the  angle  between  Ns,  the  normal  to  S,  and  Nc,  the  normal  to  C. 

3.3  Elimination  Theory 

Elimination  theory  is  a  branch  of  algebra  which  deals  with  eliminating  variables  from  equations. 
It  is  especially  useful  for  determining  when  a  system  of  equations  has  a  root.  Let  us  start  with 
the  simplest  case  -  two  polynomials  in  one  variable,  p  =  +  ...  +  po,  and  q  = 

qmx”"  +  +  ...  +  go. 

To  compute  the  resultant  of  p  and  g,  one  first  constructs  an  (n  +  m)  x  (n  +  m)  matrix  as  follows. 
Its  first  row  consists  of  p’s  coefficients,  followed  by  zeros.  The  second  row  is  obtained  by  translating 
the  first  one  to  the  right,  etc.  When  this  can  be  done  no  more,  the  same  process  is  repeated  for  q^s 
coefficients.  The  resultant  is  equal  to  the  determinant  of  this  matrix.  For  instance,  the  resultant 
of  ax'*  +  bx^  +  cx^  +  dx  +  e  and  Ax^  +  Bx^  +  Cx  +  Z)  is  the  determinant  of 

^abcdeOO^ 

0  a  b  c  d  e  0 

0  0  a  b  c  d  e 

A  B  C  D  0  0  0 

0  A  B  C  D  0  0 

0  0  A  B  C  D  0 

\0  0  0  ABCD) 

A  basic  result  in  elimination  theory  is  that  the  resultant  is  equal  to  zero  if  p  and  q  have  a 
common  root.  / 

It  is  also  possible  to  eliminate  variables  from  polynomial  systems  with  more  equations  than 
variables.  For  example,  if  we  have  three  polynomial  equations  with  two  variables,  there  is  an 
expression  in  the  coefficients  of  these  polynomials  which  is  zero  if  the  system  has  a  solution.  In 
general,  elimination  is  a  difficult  problem,  and  it  is  not  always  possible  to  explicitly  write  down 
these  expressions. 

4  The  General  Method 

In  this  section,  a  general  overview  of  our  method  for  deriving  differential  and  semi-differential 
invariants  for  curves  lying  on  surfaces  is  provided. 

We  wish  to  find  conditions  on  the  curvature  and  torsion  of  a  curve  C  which  will  allow  us  to 
determine  if  it  could  possibly  lie  on  a  certain  geometric  object  OBJ,  which  is  described  by  a  generic 
implicit  equation,  P(x,y,z)  =  0. 
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The  method  by  which  these  conditions  is  derived  proceeds  as  follows.  First,  we  use  the  local 
canonical  form  to  write  down  an  expression  for  C  in  the  vicinity  of  a  point  M  we  have  measured 
on  OBJ;  we  also  assume  that  we  have  measured  k,  t,  and  their  derivatives,  as  well  as  the  Frenet 
trihedron  at  M.  These  are  all  determined  from  the  derivatives  of  C;  so,  if  we  have  accurate 
measurements  for  C  in  the  vicinity  of  M,  we  can  directly  calculate  them.  Since  k  and  r  do  not 
depend  on  the  pose  of  the  C,  we  are  allowed  to  translate  and  rotate  OBJ  and  the  curve  on  it,  thus 
obtaining  a  new  curve  C.  Denote  the  rotated  and  translated  object  by  OBJ  new 

Every  condition  on  r  and  k  that  we  derive  for  C  is,  of  course,  also  a  condition  for  C.  The  reason 
we  apply  rigid  transformations  to  OBJ  is  because  these  allow  us  to  make  assumptions  about  C’s 
Frenet  trihedron  which  result  in  simpler  calculations;  this  wiU.  be  explained  in  the  sequel.  Let 
V{x,y,z)  be  the  implicit  equation  defining  OBJ  new 

Next,  we  substitute  C’s  local  canonical  form  into  'P{x,y,z);  This  results  in  a  Taylor  series 
in  $.  This  series  has  to  be  identically  zero,  because  C  is  contained  in  OBJ  new,  and?  therefore, 
has  to  satisfy  the  equation  which  defines  OBJ  new  This  gives  us  a  set  of  equations,  one  for  each 
coefficient  in  the  Taylor  series.  Next,  we  eliminate  from  these  equations  everything  but  C’s  curvature 
and  torsion.  For  one  curve,  we  usually  have  to  eliminate  the  Frenet  trihedron.  For  two  curves, 
we  wHl  show  that  the  Frenet  trihedrons  are  known  and  therefore  need  not  be  eliminated.  In  both 
cases,  the  elimination  gives  an  expression  that  has  to  be  zero;  and  this  is  the  sought  invariant. 

We  now  proceed  to  apply  this  paradigm  to  specific  objects. 


5  The  Case  of  a  Sphere 


In  order  to  derive  a  differential  invariant  for  a  curve  c(s)  to  lie  on  a  sphere,  we  need  to  use  only  the 
following  part  of  c’s  local  canonical  form: 


c{s)  =  c(0)  +  {s-  —)T  +  (—  +  —)N  -  —ktB  +  o( 


(5) 


Since  translation  and  rotation  do  not  change  the  curvature  and  torsion,  we  may  assume,  without 
loss  of  generality,  that  the  point  M  at  which  our  measurements  of  k  and  r  were  taken  is  at  the 
origin,  and  that  the  sphere  lies  on  the  XY  plane.  Hence  the  sphere’s  equation  is 


-\-[z  -  RY  -  =  0 


(6) 


Let  us  also  assume,  without  loss  of  generality,  that  the  sphere  has  been  rotated  so  that  T  =  (1, 0, 0) 
(see  Figure  4). 

Since  iV  is  a  unit  vector  perpendicular  to  T,  it  has  to  be  of  the  form  N  =  (0,  cos  a,  sin  a)  for 
some  a;  also,  S  =  T  X  iV  =  (0,  —  sin  a,  cos  a). 

Note  that  the  rigid  transformation  applied  to  the  sphere  has  reduced  the  Frenet  trihedron  to 
a  trihedron  depending  only  on  the  single  parameter  a.  This  is  important,  because  we  have  to 
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X  M  T  X 

Figure  4:  Rotating  and  translating  the  sphere. 


eliminate  the  trihedron,  in  order  to  obtain  a  condition  depending  only  on  n  and  r;  and,  in  general, 
the  more  variables  we  have  to  eliminate,  the  more  equations  are  necessary,  and  there  is  a  danger 
that  the  solution  will  be  extremely  complicated. 

Substituting  these  T,  JV,  B  in  Equation  5  gives  the  following  expressions  for  the  components  of 
c{s): 


,  ,  s^k  cos  a  s^(k  cos  a  +  kr  sin  a) 

!/W  =  — 2— + - 5 - 

.  .  sin  a  sin  a  —  kr  cos  a) 

=  +  - 6 - ^ 

Substituting  these  expressions  into  the  equation  of  the  sphere  (6)  gives  a  Taylor  series  in  s, 
which  has  to  be  identically  zero;  therefore  all  its  coefficients  are  zero.  The  expression  is  rather 


complicated,  so  we  don’t  write  it  down  here;  However,  its  constant  and  linear  coefficients  are 
identically  zero.  Its  coefficient  of  is 


1  -  A:(sin  a)R  =  0 

Its  coefficient  of  is 

K  (sin  a)R  —  kr  (cos  a)R  =  0 
From  (7)  and  (8),  together  with 

sin^  a  +  cos^  a  —  1  =  0 
11 


(7) 

(8) 

(9) 


(viewing  them  as  algebraic  equations  by  treating  sin  a  and  cos  a  as  algebraic  variables),  we  can 
eliminate  sin  a  and  cos  a,  to  obtain  the  identity 


B?  = 


+  {nf 


(10) 


This  gives  us  a  differential  invariant  for  a  curve  lying  on  a  sphere;  namely,  the  expression 


+  (n) 


has  to  be  a  constant.  Note  that  we  can  immediately  extract  the  sphere’s  radius. 

It  should  be  noted  that  this  condition  has  been  derived  before,  using  other  methods  (see,  for 
instance,  [6],  page  25).  We  have  nonetheless  decided  that  it  would  be  worthwhile  to  show  how  it 
is  derived  by  using  the  local  canonical  form  and  elimination  theory.  This  derivation  will  hopefully 
make  it  easier  to  follow  the  derivation  of  differential  invariants  for  curves  on  the  cylinder  and  cone, 
presented  in  the  following  sections. 


6  The  Case  of  a  Cylinder 
6.1  One  Curve,  Known  Radius 

We  now  proceed  to  derive  differential  invariants  for  a  curve  which  lies  on  a  cylinder.  To  the  best 
of  our  knowledge,  such  invariants  have  not  been  derived  before.  The  method  is  similar  to  the  one 
used  for  the  sphere;  however,  the  mathematical  details  are  considerably  more  complicated. 

Given  a  point  Af  on  a  curve  which  lies  on  a  cylinder,  we  can  assume  without  loss  of  generality 
that  the  cylinder  has  been  translated  and  rotated  so  that  M  is  at  the  origin,  and  the  cylinder  lies 
on  the  XY  plane  (recall  that  this  does  not  alter  the  curvature  and  torsion).  Let  us  further  assume 
that  it  has  been  rotated  by  some  angle  /?  so  that  the  tangent  vector  at  M  is  aligned  with  the  X-axis 
(see  Figure  5). 

Hence  T  =  (1,0,0),  and  the  cylinder’s  equation  becomes 

(x  cos/3  +  ysm/3)^  +  {z-  Rf  -  =  0  (11) 

As  in  the  case  of  the  sphere,  it  follows  that  N  =  (0,  cos  a,  sin  a)  for  some  a,  and  B  = 
(0,- sin  a,  cos  a). 

We  now  substitute  T^N into  the  local  canonical  form  (2).  This  gives  the  following  expressions 
for  the  components  of  c(s): 

^  ^  6  8^  120 


12 


Figure  5:  Rotating  and  translating  the  cylinder. 


.  .  s^K  cos  a  (k  cos  a KT  sin  a) 

y(‘)  =  - e - ^  + 

s^{k  cos  a  +  2K  r  sin  a  —  cos  a  +  kt  sin  a  —  cos  a) 

24 

o5 

S  ,  w  n  .  «  '  '  .  '  2 

cosa  +  3/c  rsmo  +  SKr  sina  -  3k  r  cosa - 
DK  K  cosa  —  KT  Sin  a  +  KT  sina  — Sktt  cosa  — kt  sma) 


.  ,  since  5®(k' sin  a  —  KT  cos  a) 

-w=^— +-i - ^ 

5^(k  sin  a  -  2  k  r  cos  a  —  sin  a  —  kt  cos  a  —  «r^  sin  a) 


(k  sin  a  —  3  k  r  cos  a  —  ok  r  cos  a  —  3  k  r  sin  a  — 

120^ 

2^»  3  3\ 

6  K  K  Sin  a  +  K  r  cos  a  —  kt  cos  a  -  3  kt  r  sin  a  +  kt  cos  a) 

Substituting  these  into  the  cylinder’s  equation  (11)  gives,  as  before,  a  Taylor  series  in  s  which 
has  to  be  identically  zero.  This  expression  is  huge  and  we  do  not  write  it  down  here;  we  need  only 
the  coeflSicients  of  the  powers  of  s  between  0  and  5. 

The  coefficients  of  the  constant  and  linear  terms  are  identically  zero. 

For  the  other  terms,  we  obtain  the  following  expressions,  after  substituting  cos  a  =  a:,  sin  a  = 
y,  cos  /?  =  2:,  sin  /?  =  w: 


For  the  coefficient  of  s'^ 
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For  the  coelRcient  of 


—  2  nyR  +  2z^  =  0 


(12) 


6zwkx  -2k  yR+2KTxR  =  0  (13) 

For  the  coefficient  of  s'* 

8 zwKT y  -2k"  yR  +  4K  txR  +  2  K^yR  +  2kt'  xR  +  2  KT^yR  -  (14) 

Q  K^X^Z^  +  8  ZWK  X  —  8  Z^K^  +  6  K^  =  0 

In  the  sequel,  it  will  be  beneficial  to  use  a  simplified  version  of  (14).  Note  that  we  can  subtract 
from  (14)  the  product  of  (12)  by  an  appropriate  constant,  and  eliminate  the  coefficient  of  in  (14) 
(it  already  has  a  yR  term,  so  we  are  not  adding  anything).  Similarly,  we  can  subtract  from  the 
new  equation  an  appropriate  multiple  of  (13),  to  remove  from  it  the  term  with  the  monomial  xzw 
—  also,  without  adding  anything  new,  as  the  set  of  monomials  of  (14)  contains  that  of  (13).  After 
grouping,  we  can  write  the  simplified  (14)  as 

Aq  +  Aix^z"^  +  A2XR .+  AzyR  +  A^yzw  =  0  (15) 

Note  that  we  can  easily  compute  the  A,’’s  as  functions  of  k  and  t.  Hence  (15)  is  equivalent  to 
(14),  but  much  simpler.  This  will  turn  out  to  be  useful. 

For  the  coefficient  of  s®,  we  obtain  the  equation 

-2  k"  yR  -  20  K?yT  xz^  -  20  kk  x^z^  -  30  zwk^x  -  10  zwkt^x  +  10  zwkt  y  + 

20 zwK  r  t/  +  10  zwk"  x  -2k^txR  +  6 k"  txR-80  z^kk  +  6  k'  v'^yR  +  (16) 

VIk^k  yRA&K  T  xR-\-2kt"  xR-2kt^xR  +  2Qkk  -\-8ktt'  yR  =  0 

in  addition  we  have 

-  1  =  0  (17) 

-  1  =  0  (18) 

Assume  now  that  the  radius  R  is  known.  In  that  case,  we  have  to  eliminate  Xfy,z,w  from 
Equations  (12,13,14,17,18)  (note  that  we  need  at  least  five  equations  in  order  to  eliminate  four 
unknowns).  All  our  attempts  to  directly  do  this,  using  various  packages  for  symbolic  computation, 
have  failed,  however,  it  is  possible  to  proceed  as  follows.  First  solve  the  system  consisting  of  the 
four  simplest  equations  (12,13,17,18);  then  substitute  the  solution  into  (14). 
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Using  the  Maple  symbolic  computation  program,  it  was  possible  to  find  a  solution  for  Equations 
(12,13,17,18).  This  solution  uses  an  auxiliary  polynomial  that  we  denote  by  pi(^): 


Pi{6)  =  81  +  (18  R?k^k^  —  18  +  81  k®  —  162  k^R?)S^  + 

36  k^RS^k  t  +  (—81 K®  +  81  k^R^  -  36  R^k'^k^  +  k^t^R^  +  2  + 

R?k^  +  18  k^t^R^)8^  -  36  k^RS^k  r  +  (18  R^k^k^  —  2  k^t^R'^k^  —  2  R?k*)6^  + 

Let  q  be  a  root  of  pi(^).  Then  the  solution  of  (12,13,17,18)  is 


z  =  R 


i 


9  —  9 

q{9qK^  +  2n  tR) 


—9 K  Rn^q"^  +  k  Rq^^n^r^  +  Rip-  -  9n  Rn^q^  —  k'^R  +  9 q 


w  = 


4^  .,3 


3,^(9  +  2  k'  r 


tR) 


y  =  - 


R{9  kY  -  ^  -  9  K^q'^  -  K^) 

g(9 QK^  +  2k  T R)k 


X  =  q 


Substituting  these  expressions  into  (14)  and  simplifying,  we  obtain  the  following  identity: 


(18  k^T^R^  —  45  R^k  +  27  R^k  k^  +  162  —  81  k^)^  +  (27  k^r  R  —  hArn  Rk^)q^  + 

(—36  k^T^R?  —  k^T^R?K^  +  6  r  k^r  +  3  k'^R^k"  k  —  2  k^r'^R?  —  5  R}k*  +  72  R^n^k*  + 
162  A;®  -  3  t^R?k  k^  -  162  k^R?  -  27  R^k  A;®)?'*  +  90  r  /c'  Rk^q^  +  (-3  k^R?k  A:  +  6  R^k"^  + 


9k^T^  R?  —  27  R?K^k‘^)<^  —  R?k^  =  0 


Let  us  denote  this  polynomial  by  ^2(9)- 

Now  we  know  that  pi{)  and  P2Q  Tnust  have  a  common  root]  therefore,  their  resultant  must  be 
zero.  This  resultant  is,  therefore,  an  invariant  for  a  curve  lying  on  a  cylinder. 

Recalling  the  definition  of  the  resultant  of  two  polynomials  (Section  3.3),  we  can  write  down 
the  resultant  of  pi()  and  P2()-  It  is  a  determinant  whose  elements  depend  on  the  curvature  and 
torsion;  if  the  curve  lies  on  a  cyhnder,  this  determinant  has  to  be  zero,  and  this  is  an  invariant  for 
a  curve  lying  on  a  cylinder. 
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Figure  6:  Logarithm  of  the  resultant  of  pi()  and  P2{)^  as  a  function  of  i?,  for  a  point  sampled  from 
a  curve  on  a  cylinder  of  radius  2. 


6.2  One  Curve,  Unknown  Radius 

6.2.1  Search  Numerically  for  the  Correct  Radius 

Suppose  we  do  not  know  the  radii  of  the  cylinders  in  the  data-base. 

There  are  two  ways  to  proceed.  We  can  simply  use  the  trivial  observation  that,  if  we  substitute  the 
correct  R  into  piQ  and  P2()?  we  will  get  two  polynomials  whose  resultant  is  zero.  We  can  therefore 
conduct  a  simple,  one-dimensional  search  for  R  which  minimizes  this  resultant. 

Experience  has  shown  us  that  this  simple  numerical  algorithm  works  quite  well.  For  example, 
Figure  6  shows,  a  plot  of  the  logarithm  of  the  resultant,  for  values  of  k  and  r  measured  on  a  curve 
on  a  cylinder  with  a  radius  of  2.  We  can  clearly  see  a  strong  minimum  at  the  correct  radius. 

6.2.2  Solve  for  the  Correct  Radius 

The  second  method  for  the  case  in  which  the  radius  is  unknown  is  to  eliminate  R^x^y^z^w  from 
Equations  (12,13,15,16,17,18).  This  can  be  done  by  solving  Equations  (12,13,15,17,18)  and  substi¬ 
tuting  the  solution  in  Equation  (16);  if  this  gives  zero,  it  means  that  these  six  equations  have  a 
common  solution,  which  is  a  necessary  condition  for  the  curve  to  lie  on  a  cylinder.  This  is  why  it 
was  important  to  define  Equation  (15),  the  simplified  version  of  (14);  we  could  not  find  a  reason¬ 
able  solution  with  (14).  However,  it  turns  out  that  Equations  (12,13,15,17,18)  do  have  a  relatively 
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simple  solution,  expressed  as  follows,  x  is  the  root  of  the  following  equation: 


81  +  {—162  —  54rA^  Aj  kJ  +  162 k^Ai  Aq  +  162^5  Ai  + 

(108  Ai  kJ  —  54  k^Ao  r  +  81  As^/c®  +  162kJ Ao  Aj  —  18  Ao  t^k^Ai  +  81  k^Aq"^  — 

54  k  A^  As  k®  —  54  rA^  As  k®  -  324  Aj  Ai  /c^  +  Or^A^^/c®  +  9  k"^  A^'^k^  —  324k®A2  Ao  + 

81  Aj^k^  +  18  k^^Aj  k  Ao  +  81  Aa^/c®)^®  +  (162 «:®A j  Ao  —  18  Ao  r^/c® Aa  + 

18k^Ao^k'^  +  108  rA^  Aa  k®  -  162  Aa^/c®  —  ISr^A^^/c®  +  18  Ao  t^k^Ai  +  6k;®A^  tk^  Aq  — 
36k^Ai/c  Ao  —  18  Ao^r^K®  +  108/c^Ao  A^  r  —  36k'*k  r  Ao  As  +  18K®Aa  k^^Ao  + 

108  K  A^  As  K®  —  324k^Ao  Aa  +  6  Ao  t^k^A^  +  162  Aa  Ai  k^  —  27 k  ^A^^k^  —  54  r  A^  Aj  — 

81  As^/c®  —  162  /c®Ao^)i5®  +  (36  k'^/c  r  Ao  Ag  —  54  k  A^  Aa  /c®  —  12  «:®A^  t  k  Ao  + 

27 +  18  Ao^r^/c®  +  Ao^t^k"^  +  Ot^A^^k®  —  54/c^Ao  A.^  r +  /c'^Ao^  +  81  k^Aq^  + 

162  kJ Ao  As  +  81  Aa^K®  +  2  Ao^t^k^k  —  54  rA^  Aa  k®  —  6  Ao  r®K®A^  — 

36  K®Aa  Ao  +  18  Ao  7-^k;®  Aa  -  36  k^Aq^k^  +  18  /c‘*A  j  /c’^Ao)^'*  +  (— 2  Ao^r^fe^/c*^  + 


/2 


/2 


;2 


/2 


6k:  A^tk  Ao+18k  Ao  k  +18k  AaK  Ao— 9k  A^  k"* — 
2  k^^Ao^)5^  +  k'^Ao^ 


Note  that  this  is  really  a  sixth-degree  equation,  as  oidy  even  powers  of  6  appear.  After  x  is  solved 
for,  we  can  easily  extract  y  from  Equation  (17).  Then,  after  substituting  the  known  values  of  x 
and  y  in  Equations  (12,13,18)  we  can  solve  for  the  remaining  unknowns  w,  z  and  R: 

Let 

e  =  K^T^x^  —  2kT  xyK  +  y^K^  +  9  K^x^y^ 

Then 


,  Root{e6'^ -9K^x‘^y‘^){kTx -yn)  „  r2  a  4  2  2^  r,  9  K^yx"^  ^ 
{it)  = - - SK^xy  - Root{e6^  -9  Krx^y^),R- j — } 


(19) 


(by  Root  of  an  equation,  we  mean  the  root  of  the  equation  when  viewed  as  an  equation  in  6).  The 
equations  in  (19)  are  trivial  to  solve  by  taking  square  roots. 

The  reader  may  ask  why  we  did  not  apply  this  trick  to  simplify  the  solution  of  the  equations 
for  the  case  in  which  the  radius  is  known.  If  we  had  done  that,  it  would  not  have  been  possible  to 
obtain  a  function  of  k  and  r  alone;  the  x  would  stiU  have  been  there!  And,  as  long  as  it  is  there, 
we  cannot  find  a  condition  on  k  and  r,  as  desired,  but  only  a  condition  on  x,  n  and  r. 

There  is  also  a  direct  solution  to  the  system  of  equations  (12,13,15,17,18),  in  which  all  the 
unknowns  R,  x,  y,  z,  w  are  written  in  terms  of  k,  r,  and  their  derivatives;  however,  that  expression 
is  truly  enormous,  covering  three  fuU  pages  even  in  a  small  format!  For  all  practical  purposes, 
it  is  better  to  use  the  solution  above,  which  first  extracts  x  and  y  and  then  solves  for  the  other 
unknowns. 
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6.2.3  Comparison  of  Methods 

While  the  second  method  is  straightforward  and  does  not  require  any  search  (as  opposed  to  the 
first  method),  it  has  the  drawback  of  requiring  the  fifth  derivative  of  the  curve,  which  appears  in 
Equation  (16).  Calculating  the  third  derivative  of  the  curvature  and  the  second  derivative  of  the 
torsion  requires  the  fifth  derivative  of  the  curve.)  The  first  method  requires  a  numerical  search  for 
the  correct  radius,  but  uses  only  the  fourth  derivative  of  the  curve.  Depending  on  how  accurate 
the  measurements  are,  one  may  opt  for  using  the  first  or  the  second  method. 

6.3  The  Case  of  a  Cylinder  with  Two  Intersecting  Curves  and  Unknown  Radius 

Suppose  we  have  two  curves  on  the  cylinder,  intersecting  at  a  point  M.  For  instance,  one  can  design 
a  sensory  robot  to  traverse  a  point  twice,  in  two  different  directions.  Another  possible  source  is  an 
intersecting  pattern  of  structured  light  rays.  It  turns  out  that  a  particularly  simple  invariant  can 
be  written  in  this  case. 

We  refer  to  the  curves  as  “first”  and  “second”  (it  makes  no  difference  which  is  which,  of  course). 
As  noted  in  Section  3.2,  two  intersecting  curves  on  the  surface  allow  us  to  compute  its  normal  N 
(denoted  this  way  to  prevent  confusion  with  N,  the  normal  to  a  curve).  We  may,  as  before,  translate 
and  rotate  the  cylinder  so  that  the  intersection  point  M  is  at  the  origin,  the  cylinder  lies  on  the  XY 
plane,  and  the  tangent  vector  to  the  first  curve  is  (1,0,0).  The  difference  is  that  now,  as  opposed  to 
when  we  only  had  a  single  curve,  we  know  the  normal  N  and  the  binormal  B  of  the  new  curves;  this 
is  because  we  now  know  that  the  rotation  and  translation  not  only  move  M  to  the  origin  and  align 
the  tangent  to  the  first  curve  with  (1,0, 0),  but  also  align  Af  with  (0, 0, 1).  Let  us  look  at  the  triplet 
{T,N ,B)  for  the  first  curve  (before  the  rotation).  We  can  calculate  the  inner  products  {N,T)  and 
{N,Af).  These  inner  products  do  not  change  after  the  rotation  of  the  curve;  if  N  is  rotated  into 
Nnew,  then,  since  T  is  rotated  into  (1,0,0),  we  have  the  equality  (iV„e«,,  (1,0,0))  =  {N,T),  and, 
since  Af  is  rotated  into  (0,0,1),  we  have  the  equality  (iV„ev;,  (0,0, 1))  =  {N,Af).  Since  is  a 
unit  vector,  we  can  recover  it;  and,  since  we  know  the  tangent  and  normal  of  the  new  curve,  we 
know  its  binormal,  which  is  equal  to  their  vector  product.  Following  a  similar  argument,  we  also 
know  the  Frenet  trihedron  of  the  (new)  second  curve  . 

As  before,  let  /?  denote  the  angle  at  which  the  cylinder  is  aligned  relative  to  the  XY  plane. 
Let  us  denote  the  tangent,  normal  and  binormal  of  the  first  curve  by  (1,0,0),  (0,cosQ,sino)  and 
(0,  -  sin  a,  cos  a),  and  those  of  the  second  curve  by  {Ti,T2, 0),  (JVi,  N2,  W3),  {Bi,B2,  B3)  (remember 
that  aU  these  coordinates  are  now  known).  Note  that  the  ^-coordinate  of  both  tangents  has  to  be 
zero,  as  they  are  both  in  the  tangent  plane  which,  after  the  rigid  transformation,  is  the  XY  plane. 

Substituting  these  expressions  into  the  local  canonical  form,  then  into  the  cylinder’s  equation, 
and  equating  coefficients  to  zero,  results  in  the  following  equations  (A  is  the  curvature  of  the  second 
curve): 
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For  the  coefficient  of  5^,  first  curve,  we  have 


2  2^  —  2  A:  sin(a)^  =  0  (20) 

For  the  coefficient  of  s^,  first  curve,  we  have 

Q  zwk  cos(a)  —  2  K  sm{a)R  + 2  kr  cos{a)R  =  0  (21) 

For  the  coefficient  of  s^,  second  curve,  we  have 

2  -  2  KN3  R  +  2  z'^Ti^  -  2  T^^z"^  +  4  2T1  WT2  =  0  (22) 

Also 

z^  —  1  —  0  (23) 

(where,  as  before,  cos/3  =  z,sin/?  =  w). 

Eliminating  w^z  and  R  from  (20,21,22,23)  results  in  the  identity 


9  Ti^K^sin^cccos^a  —  6  r  T2  TiK^cos^a  sin  a  —  9  KN^K^cos^a  sin  a  +  K^r^Ta^cos^a  + 

6  k^Ti  T2  K  sin^o:  cos  a  —  2  kk  t  T2^  sin  a  cos  a  +  =  0 

(remember  that  a  is  known,  and  does  not  have  to  be  eliminated). 

This  is  an  invariant  for  two  intersecting  curves  which  can  be  used  to  test  whether  they  lie  on 
a  cylinder.  The  invariant  depends  on  the  curvature  and  torsion  of  one  curve,  and  the  curvature  of 
the  other;  therefore,  it  does  not  require  any  derivatives  of  order  higher  than  3. 

7  The  Case  of  a  Cone  with  Two  Intersecting  Curves 

We  have  not  addressed  the  problem  of  finding  invariants  for  a  cone  using  a  single  curve;  because 
a  cone  has  more  degrees  of  freedom  than  a  sphere  or  a  cylinder,  this  would  necessitate  using  the 
sixth  derivative  of  the  curve  to  express  such  an  invariant. 

We  win  show  in  this  section  how  two  intersecting  curves  yield  an  invariant  for  the  cone.  We  will 
not  go  into  aU  the  details,  as  the  method  resembles  the  one  used  for  a  cylinder  with  two  intersecting 
curves. 

First,  the  cone  is  rotated  and  translated  so  that  its  apex  is  at  the  origin,  and  the  point  of 
intersection  of  the  two  curves,  M,  lies  on  the  XY  plane,  which  is  also  the  tangent  plane  at  M. 
Then,  it  is  rotated  in  the  XY  plane  so  that  the  tangent  vector  of  the  first  curve  is  (1,0,0).  As  we 
did  for  the  cyfinder,  we  can  extract  the  tangent,  normal,  and  binormal  vectors  to  the  two  curves  at 
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Figure  7:  Rotating  and  translating  the  cone. 


their  new  locations;  denote  the  normal  to  the  first  curve  at  M  by  (0,cos/3,sin^).  Note  that  now 
M  does  not  lie  at  the  origin,  but  at  an  (unknown)  distance  j/q  from  it.  The  (unknown  -  as  for  the 
cylinder)  rotation  angle  of  the  cone  in  the  XY  plane  is  denoted  by  a,  and  the  (unknown)  rotation 
angle  around  the  Y  axis  (Figure  7)  is  denoted  this  is  just  half  of  the  cone’s  opening  angle. 

It  is  then  a  trivial  matter  to  write  down  the  equation  of  the  rotated  and  translated  cone,  and  to 
substitute  into  it  the  local  canonical  forms  of  the  two  curves.  As  before,  the  coefficients  of  the  two 
resulting  Taylor  series  have  to  be  zero,  resulting  in  the  following  equations.  Let  K  and  T  denote 
the  curvature  and  torsion  of  the  second  curve,  and  Ti,  etc.  the  components  of  its  Frenet  trihedron 
vectors.  S  stands  for  sin/?,  C  for  cos^,  x  for  cos^,  y  for  sinO,  z  for  cos  a,  w  for  sin  a  (remember 
that  S  and  C  are  known,  and  do  not  have  to  be  solved  for). 

For  the  coefficient  of  5^,  first  curve,  we  have 


2x^z^  +  2xkSyoy  =  0 
For  the  coefficient  of  5^,  first  curve,  we  have 

—  2xkTCyoy  +  6yxwkS  +  2xk  Syoy  —  Gx'^zwkC  =  0 

which  can  be  written  more  compactly  as 

Ai  x'^zw  +  A2  xyyo  +  A3  xyw  =  0 
For  the  coefficient  of  5^,  second  curve,  we  have 

-4x^zwTiT2  +  2x'^T2‘^ -2x^T2'^z'^ +  2xKN3yoy  =  0 


(24) 


(25) 


(26) 


(27) 
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Just  as  for  the  cylinder,  we  can  subtract  from  (27)  appropriate  multiples  of  (24)  and  (25),  and 
obtain  the  simpler  form 


A4  +  As  xypo  +  ^6  =  0  (28) 

For  the  coefficient  of  s^,  second  curve,  we  have 

6  yxzKNz  T2  +  6  yxwKNz  Ti  +  6  x^Tz  KN2  -  6  x^T^  KNzz^  —  2  xKT  BzyoVA 
2xKlNzyoy-Qx^zwTiKN2-&x'^zwKNiT2A&x'^z^TiKNi  =  0  (29) 

which,  as  before,  can  be  reduced  to 

At  xyw  +  As  xyz  +  Ag  +  Aio  xyyo  =  0  (30) 

Note  that  the  A,’s  can  be  readily  computed  from  the  known  quantities  -  the  curvature,  torsion, 
and  Frenet  trihedron  of  the  two  curves. 

We  also  have  the  equations 


x^  +  y^-1  =  0  (31) 

and 

w^  +  z^  --1  =  0  (32) 

It  is  possible  to  eliminate  {x,  y,  z,  yo}  from  these  six  equations  (24,26,28,30,31,32),  and  obtain 

—2  SkA4^ As  Aio  A2  A7  —  A^^ A2^ Ar^ A^  +  2  SkA4^ A2  A7  As  As  Aq  —  SkAq  A^^ A2^A7  As  ~ 

As^ A^ A^^ As  —  +  SkA^  A^^ A^^ As’^  —  As^ A4  A^ A^  ylio  — 

S’^k’^As^As  A9  >l2  ^8  +  2  S'^k^As  A4^A2  As  As  A9  -  S^k^A4^Ai^As^A5  - 

S^K^Ae  A4^Ai^A7  As -2  S\^A4^Ai  A5  As  ^3  ^9  +  2  S^k^Aq  A4^Ai  A7  As  A9  + 

SkAq  A4^ As  ^10  ^2  -4-8  —  2  iS/cu44^^3^-Aio  -A5  -A9  +  kP A4^ Ai^ As^  H-  Aq  A4^ Ai  As^ A2  + 

S^k^A4^Ai^A7^  +  SkA4^A2^As^  +  S^K^Ae^A4  As^Ag^  +  2  A4  As^Aio  - 

2  A4  As  Ag  As^ A2  A7  —  3  SkAq  A4  As  Ag  As  A2  As  +  S/Cj44^A3^Aio^  +  2  SKA4  As  Ag  As^ Ai  As  + 

2  A4^ Ai  As  As  Aiq  —  2  SkA4^ Ai  As  As  Aio  ^5  +  Aq^  A4  As  Ag  A2  A7  +  SkA4^ A2^ A7^  — 

Aq^ A4  As  Ag  Ai  As  —  SkAq^ As  Ag  A2  A7  As  A  SkAq  As  Ag  Ai  As^ A7  H- 
SkAq^As^  Ag  Aio  As  A  Aq^  As  Ag  Ai  As  As  —  2  A4^  Ai^ A7^ As  + 

Aq  A4^ Ai  A7^ A2  ~  Aq  A4^i4i  A7  As  Aio  ■“  -44^-A3^Aio^A5  —  3  Aq  A4  Ai  A7  ^3  A9  As  + 

2  A4^ As  Aio  A2  A7  As  —  SkAq  A4  Ai  A7^A2  As  4-  SkAq  A4  Ai  A7  As  Aio  As  ” 

Aq^ A4  Ai  A7  A2  As  A  K,"^ Aq  A4  Ai^ A7  As  As  SKA4  Ai^A7^As^  ~  0 
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This  is  an  invariant  for  two  curves  on  a  cone.  It  depends  on  the  curvature  and  torsion  of  the  two 
curves;  therefore,  it  does  not  require  any  derivatives  of  order  higher  than  three.  This  is  an  invariant 
for  two  curves  on  a  cone. 

8  Numerical  Computation  of  Derivatives 

The  algorithms  suggested  here  require  computing  the  derivatives  of  a  curve  in  3D  space.  The  prob¬ 
lem  of  computing  high-order  derivatives  from  discrete  data  was  addressed  in  [25].  The  derivatives 
at  each  point  are  calculated  by  convolving  appropriate  differentiation  filters  with  the  given  curve. 
One  way  of  deriving  such  filters  is  based  on  fitting  high-order  polynomials  to  the  data  curve  and 
differentiating  the  polynomials.  We  do  not  need  to  do  the  fitting  for  each  actual  curve;  it  is  only 
done  in  deriving  the  filters. 

In  deriving  the  filters,  the  data  curve  /  is  approximated  by  a  linear  combination  of  orthogonal 
polynomials  of  orders  0, . . . ,  Z: 

i 

Fi{-x)  -  w{x)Y^aiPi{x) 

i=0 

where  Pi{x^  are  polynomials  which  are  orthonormal  with  respect  to  a  weight  function  w{x).  The 
coefficients  a,  are  determined  by  the  condition  that  the  polynomial  fits  the  curve  in  the  sense  of 
(weighted)  least  squares.  It  can  be  proved  that  if  the  curve  /  is  a  polynomial  of  order  up  to  Z,  than 
the  above  filter  yields  an  exact  A;-th  derivatives  when  the  coefficients  a,-  are 

a,-  =  P^\q) 

In  practice,  good  results  are  obtained  for  any  reasonably  smooth  /  (not  only  polynomials),  as  long 
as  the  the  order  I  of  the  filter  is  larger  than  the  desired  order  k  of  the  derivative.  However,  a  high 
Z  requires  a  filter  with  a  wide  support. 

Discrete  versions  of  this  method  on  a  finite  interval  are  described  in  detail  in  [15].  In  particular, 
the  Krawtchouk  and  the  discrete  Chebyshev  polynomials  were  studied  and  closed  form  formulas 
for  them  were  given  up  to  fifth  order.  However,  it  was  shown  in  [25]  that  continuous  polynomials, 
defined  on  a  finite  interval,  are  just  as  effective  but  much  simpler  to  calculate.  Good  results  were 
obtained  using  the  Legendre  and  continuous  Chebyshev  polynomials. 

For  example,  see  Figure  8  for  the  derivative  of  the  spherical  invariant  (Equation  10),  for  the 
curve  plotted  in  Figure  3,  when  noise  of  variance  equal  to  5  percent  of  the  distance  between  the 
points  was  added  to  it.  Derivatives  were  computed  using  the  method  described  in  [25].  The 
derivative  is  relatively  small  for  the  part  of  the  curve  that  lies  on  the  sphere  (0  <  t  <  1),  and 
significantly  changes  when  the  curve  crosses  over  to  the  cylinder  (at  t  =  1).  Note  that  computing 
the  invariant’s  derivatives  requires  the  first  derivative  of  the  torsion  and  the  second  derivative  of 
the  curvature  —  that  is,  the  fourth  derivative  of  the  curve. 
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Figure  8:  Derivative  of  the  spherical  invariant  (Equation  10)  for  the  curve  in  Figure  3,  after  noise 
had  been  added  to  it.  Note  the  change  in  the  derivative  when  the  curve  crosses  over  from  the 
sphere  to  the  cylinder  (at  t  =  1). 


9  Semi-Differential  Invariants 

In  this  section  we  study  curve  invariants  which  use  only  curvature  (this  requires  computing  only  the 
first  and  second  derivatives  of  the  curve).  We  also  assume  that  the  only  primitives  the  recognition 
system  can  encounter  are  spheres,  cylinders,  cones,  and  tori.  When  the  information  from  one  point 
is  not  enough  to  uniquely  determine  the  object,  we  will  use  an  additional  point  or  two  on  the  curve 
to  help  disambiguate  the  object. 

9.1  Object  Recognition  from  Two  Intersecting  Curves 

Given  two  intersecting  curves  C\  and  C2,  we  extract  Ti,  Ni,  Bi,  Kci,T2.,  N2,  B2,  Kc2  at  the  inter¬ 
section  point  M.  These  are  the  Frenet  trihedrons  and  the  curvature  for  both  curves  respectively. 
Recall  that  Ns,  the  normal  to  the  surface  at  M,  equals  T\  XT2. 

For  each  curve  we  compute  0,  the  angle  between  Ns  and  the  curve’s  normal.  The  surface  normal 
curvature  equals  =  kc  cos0,  and  =  kj  sin^yS  +  «;2Cos^/3,  where  ki,K2  are  the  principal 

curvatures  of  the  surface  at  M,  and  /?  is  the  angle  between  the  tangent  to  the  curve  and  K2,  the 
second  principal  direction. 

Given  two  curves  we  have  two  equations  for  the  surface  normal  curvature  in  the  three  unknowns 
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Ki,  A;2>  and  /?; 

knsI  =  «i  sin^  /8  +  K2  cos^  /S 

«jVs2  =  «i  sm^(/3  +  <!>)  +  K2  cos^i/S  +  </>),  (33) 

where  <j>  is  the  angle  between  T\  and  T2.  Usually,  it  is  impossible  to  solve  such  a  system;  however, 
if  we  know  in  advance  that  the  geometric  primitives  can  only  be  spheres,  cylinders,  cones,  and  tori, 
it  is  possible  to  identify  them  and  extract  their  parameters. 

Sphere 

In  this  case  ki  =  K2  and  consequently  /catji  =  Kj\^g2-  For  all  other  objects  (cylinders,  cones,  and 
tori)  the  two  principal  curvatures  are  not  equal;  therefore,  the  two  normal  curvatures  are  identical 
only  in  the  degenerate  case  in  which  the  angles  between  the  curves’  tangents  and  kj  are  equal. 
Therefore,  if  the  surface  normal  curvatures  corresponding  to  the  two  curves  are  equal,  we  can 
assume  with  high  probability  that  we  are  dealing  with  a  sphere. 

The  sphere’s  radius  is  then  R  =  1/kn^,  and  its  center  is  at  Af  +  RNs-  Using  an  additional 
point  (Section  8.2),  we  can  determine  whether  the  object  is  indeed  a  sphere. 

Cylinder 

If  the  given  object  is  a  cylinder,  its  parameters  can  be  recovered  as  follows.  Since  ki  =  0,  the 
surface  normal  equations  are  reduced  to  two  equations  with  two  unknowns.  Solving  them,  we  can 
recover  K2  and  the  principal  directions  ki,K2-  The  cylinder’s  radius  is  jR  =  —  amj  tjjg  orientation 
of  its  axis  is  Kj.  A  point  on  the  axis  is 


C  =  M  +  RNs. 

It  is  important  to  note  that  this  does  not  prove  that  the  object  is  a  cylinder.  That  has  to  be 
verified  using  an  additional  point  on  the  curve  (see  Section  9.2). 

Cone 

Assume  the  object  is  a  cone.  As  for  the  cylinder,  =  0,  and  we  can  recover  K2  and  the  principal 
directions  ki,K2-  The  radius  of  the  cone  at  Af  is  i?  =  ^cosa,  where  a  is  the  cone’s  opening 
angle.  The  apex  is  located  at  Af  +  and  the  axis  orientation  is  ki  cos  a  +  k*2  sin  a.  a  can  be 

determined  from  an  additional  point  on  the  curve  (see  Section  9.2). 

Torus 

From  (33)  we  cannot  recover  the  torus,  because  the  number  of  unknowns  is  3.  We  wiU  parameterized 
our  solution  as  a  parameter  of  /?.  For  a  given  /?,  we  can  recover  ki,  K2.  The  values  of  ki,K2  change 
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on  the  torus  as  a  function  of  7,  the  angle  between  the  major  radius  of  the  torus,  J?,  and  the  vector 
to  the  current  point  on  the  torus.  Ki, ^2  are  given  as  a  function  of  7: 

-  cos  7  1 

Ki  =  — -  K2  =  - 

K  +  r  cos  7  r 

where  R  and  r  are  the  major  and  minor  radii  of  the  torus  respectively  (see  Figure  9). 


Figure  9:  Torus. 

Given  «;i,K2?7  we  can  recover  R,r  as  follows: 

R=—i - )cos7  r=  — 

Ki  K2  «2 

The  orientation  of  the  torus,  can  be  recovered  by:  Nt  =  Nssinj  +  K2  cos  7.  The  center  of  the 
torus  is  then  at 

C  —  Af  +  NtT  +  {Ns  cos  7  -  /?2  sin  ^)R 

P  and  7  can  be  determined  by  an  additional  point  on  the  curve  (see  Section  9.2). 

9.2  Verification  Using  an  Additional  Point 

The  hypothesis  about  an  object  and  its  parameters  can  be  verified  by  using  an  additional  point  on 
one  of  the  curves.  For  the  hypothesis  to  be  correct,  several  constraints  must  be  satisfied. 

•  The  point  M  must  lie  on  the  surface.  This  means  that  if  P  is  the  object’s  implicit  equation, 
P{M)  =  0. 

•  Tc,  the  curve’s  tangent,  must  be  orthogonal  to  Ns  at  the  point.  Thus 

Ns  ’Tc  =  0. 

•  If  0  is  the  angle  between  Ns  and  Nc^  then  ks  =  nc  cos  0  (see  Equation  4),  where  the  value  of 
Ks  is  determined  by  the  principal  curvatures  Ki  and  «2  ^ind  the  angle  between  them  and  Tc- 
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Therefore,  each  additional  point  yields  three  additional  equations  which  have  to  be  satisfied. 
These  equations  can  be  used  to  verify  hypotheses  or  to  determine  the  value  of  unknown  parameters. 

If  the  additional  points  are  not  on  a  curve,  and  we  don’t  have  any  differential  properties  asso¬ 
ciated  with  them,  we  still  have  the  first  condition  (they  have  to  satisfy  the  surface  equation).  In 
that  case,  we  will  need  more  points;  this  is  a  typical  tradeoff  for  semi-differential  invariants. 

9.3  Object  Recognition  from  One  Curve 

When  two  intersecting  curves  are  given,  we  are  able  to  recover  Ns  and  thus  we  know  the  angle  0 
between  Ns  and  Nc-  When  we  are  given  only  one  curve,  9  is  an  unknown  parameter  which  has  to 
be  recovered. 

Sphere 

In  this  case  ki  =  K2,  and  consequently  =  1/R.  For  every  value  of  9,  the  surface  normal  and 

the  sphere  s  radius  are  determined  as  follows,  where  {Tq,  Nc,  Be)  are  the  Frenet  trihedron  of  the 
curve: 


Ns  =  (cos  9)Nc  +  (sin  9)Bc 
We  can  then  recover  the  center  of  the  sphere,. 


R  = 


1 

KC  COS  9 


M  +  RNs  =  M  +  -f  Be  tan  9 
i^c 


(34) 


Thus  we  have  a  family  of  possible  spheres,  parameterized  by  9. 

Given  additional  points,  we  can  proceed  as  follows:  either  substitute  them  into  the  (hypothe¬ 
sized)  sphere  s  equation,  or  (if  they  are  on  a  curve)  use  the  verification  method  described  in  Section 
9.2.  Alternatively,  given  two  points  on  a  curve,  applying  Equation(34)  to  both  of  them  results  in 
four  linear  equations  in  cos ^i,  tan ^i,  cos ^2,  and  tan  ^2-  The  solution  is  verified  by  checking  if  the 
two  angles  satisfy 


tan  9i 


y/1  —  COS^  9i 
cos  9i 


Cylinder 


In  the  case  of  the  cylinder  we  know  that  ki  =  0  and  «2  =  Given  a  point  Mi  on  the  curve,  the 
two  unknowns  are  9i  and  f3i.  When  they  are  given,  the  cylinder  is  uniquely  defined.  Note  that  /cj 
is  the  axis  of  the  cylinder,  so  it  has  to  be  the  same  for  every  point  on  the  cylinder.  ’We  will  now 
use  these  facts  to  define  R  and  kI  the  axis  of  the  cylinder  as  functions  of  9i  and  /?i  (see  Equations 
33,  4): 

cos^  Pi 


R  = 


Kci  COS  ' 
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Ns  =  cosOiNc  +  sin0iJ?c 


Ki  =  Tc  sin  /?i  +  (Tc  X  Ns)  cos  /3i 


(35) 


A  point  on  the  aods  is 

Cl  =  Ml  +  RNsi 

Given  an  additional  point,  its  /?2  and  02  can  be  recovered  as  follows: 

^  /C0S^/?2x 

P2  =  arcsin(rc2  *  «i),  ^2  =  arccos( - — 

From  them  we  can  recover  the  point  on  the  cixis  C2  closest  to  the  second  point,  and  both  points 
must  lie  on  the  cylinder’s  cixis,  which  is  parallel  to  kI;  therefore 


(Cl  -  C2)  X  «i  =  0, 

which  gives  us  two  equations  in  two  unknowns,  which  can  be  solved  for  ^1  and  /?i. 

These  two  points  give  the  equation  of  the  cylinder  that  passes  through  them  and  satisfies  the 
given  constraints.  In  addition,  from  (35)  j32  -  arccos((rc'2  X  A52)  *  /^),  which  gives  an  additional 
constraint  to  verify  that  this  is  indeed  a  cylinder  with  the  computed  parameters. 


Cone 

The  case  of  the  cone  is  similar  to  that  of  the  cylinder  but  slightly  more  complicated.  Given  two 
points  on  a  curve  we  would  like  to  find  the  angles  ctiid  /?2-  These  angles  parameterize 

the  local  surface  structure  of  the  two  points.  At  first  we  will  exploit  the  fact  that  the  line  from 
the  point  on  the  surface  in  the  direction  of  must  pass  through  the  tip  of  the  cone.  Thus  we 
have  a  constraint  that  the  two  such  lines  obtained  from  the  two  points  must  intersect.  The  point 
C  =  P  +  K2NS  lies  on  the  central  axis  of  the  cone.  Therefore  we  two  additional  constraints  which 
are  due  to  the  fact  that  Ci,  C2,  and  the  tip  of  the  cone  lie  on  the  same  line.  Finally,  the  angle  of 
the  cone  a  must  be  the  same  for  both  surface  points.  As  a  is  the  angle  between  and  the  axis  of 
the  cone,  we  can  write  an  additional  constraint  enforcing  the  uniqueness  of  a.  With  the  four  above 
mentioned  constraints  we  can  recover  the  values  of  the  unknown  angles  and  recover  the  shape  of 
the  object. 

As  for  the  cylinder,  these  two  points  give  the  equation  of  the  cone  that  passes  through  them 
and  satisfies  the  given  constraints.  However,  an  additional  point  is  needed  to  verify  that  this  is 
indeed  the  real  object. 
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Torus 


In  order  to  be  able  to  recover  the  seven  parameters  of  the  torus,  we  parameterize  them  by  four 
local  parameters  of  one  point.  The  parameters  are  and  7.  As  described  above  these  four 

parameters  are  enough  to  describe  the  torus.  In  order  to  recover  those  parameters  we  need  two 
additional  points  because  each  point  yields  three  constraints.  Thus  using  three  points  we  can 
recover  the  shape  of  the  torus  and  verify  that  it  is  indeed  a  torus. 

9.3.1  Experimental  Results 

The  algorithm  for  a  single  curve  has  been  tested  on  real  data  obtained  by  the  Brown  k  Sharpe 
Company  using  their  coordinate  measuring  machines  (Figure  1).  The  data  is  a  curve  measured  on 
a  cylinder.  For  each  point  on  the  curve  Tc,Nc,Bc,  and  kc  are  estimated.  Using  the  algorithm 
described  above,  the  problem  is  reduced  to  solving  for  cos  0%  and  cos  j3i ,  where  all  other  parameters 
are  expressed  as  functions  of  these  unknown  values.  The  correct  values  must  satisfy  four  equations 
and  must  satisfy  the  constraints  that  the  absolute  values  of  the  cosines  and  sines  of  the  various 
angles  must  be  less  than  1.  The  values  of  the  unknowns  are  found  using  non-linear  least  squares 
optimization  techniques.  We  used  the  Levenberg-Marquardt  procedure  of  the  MINPACK  library 
[16]. 

We  chose  200  pairs  of  points  at  random  and  ran  the  minimization  procedure  on  them  using 
several  initial  conditions  for  the  minimization.  Even  though  the  data  is  noisy,  most  pairs  of  points 
yielded  results  close  to  the  correct  shape.  The  results  were  sorted  according  to  the  least-squares 
error  (LSE)  of  the  four  equations.  We  show  the  five  cylinders  with  the  smallest  LSEs  in  Figure 
10(a).  One  of  these  results  and  the  original  data  are  shown  in  Figure  10(b).  It  is  important  to  note 
that  only  the  data  for  the  two  points  and  their  derivatives  mentioned  above  was  used  to  recover 
the  shape  of  the  cylinder.  Additional  points  can  then  be  used,  if  desired,  to  get  a  better  estimate 
of  the  shape. 

10  Conclusions 

A  novel  method  of  recognizing  some  surfaces,  given  curve(s)  on  them,  was  presented.  It  proceeds 
by  using  invariants  which  are  computed  on  curves,  but  which  supply  information  about  the  type 
of  surfaces  the  curves  can  possibly  lie  on. 

The  method  can  use  3D  curves  derived  from  stereo  and  structured  light;  it  is  particularly  useful 
when  given  the  output  of  measuring  devices  which  produce  such  curves  (for  instance,  sensory  robots 
and  coordinate  measuring  machines). 

The  main  advantage  of  the  proposed  method  compared  to  algebraic  methods  is  in  its  loced 
nature,  which  enables  it  to  segment  and  recognize  curves  (and  the  surfaces  they  lie  on)  even  if  the 
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(a)  (b) 

Figure  10:  (a)  The  five  recovered  cylinders  with  the  lowest  LSE.  (b)  The  recovered  shape  of  the 
cylinder  and  the  data  points. 


curves  lie  on  more  than  one  geometric  primitive.  Also,  it  requires  a  far  smaller  number  of  curve 
points  than  the  algebraic  method,  to  recognize  a  single  primitive. 
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